Fractional generalization of Fick's law: a microscopic approach 
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In the study of transport in inhomogeneous systems it is common to construct transport equations 
invoking the inhomogeneous Fick law. The validity of this approach requires that at least two 
ingredients be present in the system. First, finite characteristic length and time scales associated 
to the dominant transport process must exist. Secondly, the transport mechanism must satisfy a 
microscopic symmetry: global reversibility. Global reversibility is often satisfied in nature. However, 
many complex systems exhibit a lack of finite characteristic scales. In this Letter we show how to 
construct a generalization of the inhomogeneous Fick law that does not require the existence of 
characteristic scales while still satisfying global reversibility. 

PACS numbers: 05.40.Fb, 05.10.Gg, 05.60.-k 
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Over a hundred years ago, A. Fick proposed the famous 
formula that bears his name whilst studying transport 
processes in saline aqueous solutions He observed 
that a salt flux was driven by the presence of a salt den- 
sity gradient (in the absence of other forces) and hypoth- 
esized that they were linearly related to each other. In 
one dimension, Fick's relation would read: 



T F = -A 



dn 
dx 



(1) 



where n is the salt density and the minus sign makes 
sure that the flux tends to reduce the gradient. Since 
then, Fick's law has found application in countless sys- 
tems, and expressions like Eq. {1} are applied to the 
transport of many quantities: electrons and ions, neutral 
particles, energy or momentum, risk, chemical rcactants, 
or currency, just to name a few. 

Using Fick's expression in the standard continuity 
equation yields the familiar classical diffusive equation, 



dn A® 2n 
dt dx 2 



(2) 



where n represents now the density of the quantity of 
interest. The constant A is the diffusive coefficient, which 
is related to some well-defined characteristic length (I) 
and time (r) scales of the underlying dominant transport 
process (in the case of salt transport, the mean-free-path 
between collisions with the background water molecules 
and the inter-collision time): A ~ I 2 /t. By well-defined, 
we mean that I and r are both finite and much smaller 
than the system size and lifespan, respectively. 

Many physical systems are, however, inhomogeneous. 
The usual inhomogeneous extension of Fick's law is: 



-A(x 



dn 
dx 



dn 
dt 



d_ 

dx 



A(x) 



dn 
dx 



(3) 



In spite of its familiarity, this expression is not the only 
possible extension that reduces to Eq. @ when A{x) = 



A = const. In fact, Fick's law is a particular case of the 
more general Fokker-Planck diffusive law Q , that defines 
the flux by, 



r FP = B{x)n - |- [A{x)n] 



(4) 



Eq. ([3]) requires that B{x) = dA(x)/dx. Traditionally, 
the validity of this constraint (and thus Eq. has been 
justified by invoking the presence of a particular symme- 
try of the underlying microscopic transport mechanism: 
local reversibility (LR) Q. The symmetry applies when 
the probability (per unit time) of a particle being trans- 
ported from any location x to any other location x' equals 
that of being transported from x' to x [For instance, con- 
sider reactants moving in a reactive medium. LR holds 
whenever the probability rate of moving between any x 
and x' without reacting depends on the total number of 
reactive centers found along the reactant path. Similar 
situations are very common in nature.]. The condition 
for the validity of Fick's law is however less stringent 
than LR: global reversibility (GR) is sufficient. Namely, 
that the probability (per unit time) of a particle leaving 
x be the same as the probability (per unit time) of arriv- 
ing at x. Clearly, LR implies GR but the converse is not 
true 0). 

The purpose of this Letter is to revisit this discussion 
in a different but related context, which has emerged as 
an important paradigm in the last decade: that of scale- 
free (or fractional) diffusion 0, Q . Scale- free transport 
appears in systems in which some of the transport char- 
acteristic scales, on which Fick's law is fundamentally 
based, are absent. The question we will try to answer 
in this Letter is: which is the expression that correctly 
represents fluxes in an inhomogeneous system satisfy- 
ing GR but lacking such characteristic scales? The rel- 
evance and timeliness of the answer is justified by the 
recent surge of interest in physical, economical, biolog- 
ical and social systems in which scale-free transport is 
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observed [1, 0]. An example is provided by those sys- 
tems in which (particle or energy) transport takes place 
via correlated avalanches. This seems to be the case of 
some magnetically confined fusion plasmas [EJ1, 0] , the 
propagation of forest fires @], earthquakes [9|] or solar 
flares [10(. Transport events in these cases have a maxi- 
mum size that is only limited by the system size L, and 
therefore a characteristic size that diverges with (some 
power of) L. In some of these systems, LR is satisfied at 
least in certain limits. For instance, in the propagation 
of forest fires, the transition rate at which fire propagates 
from tree A to tree B would heavily depend on how many 
active trees exist between them, thus satisfying LR. 

The discussion will proceed in parallel to the previ- 
ous one on classical diffusion. We start by reviewing the 
scale- free version of the homogeneous Fick law 0, 0| • The 
relevant transport equation is then expressed in terms 
of fractional differential operators (FDOs) [ll[. To un- 
derstand why, we consider the simplest one, known as 
the Markovian, symmetric fractional diffusion equation 
(MsFDE): 

dn d a n . 

m =Aa W¥' ae(0 ' 2] ' (5) 

where A a is a constant. The Riesz FDO is defined as 
d a -1 



d\x\ a 2cos(?ra/2) 



(6) 



where the Ricmann-Liouvillc FDOs of order a are [ll| : 

i d m r f{x') 



aD a x f 
a D a J 



F(m — a) dx r 
(-1)"' d m 
T(m — a) dx r 



(x 



r f\a—m-\- 



/(*') 



\a—m-\- 



T dx', (7) 



being m the integer number verifying that m — 1 < 
a < m. In spite of their complicated appearance, the 
theory of fractional differential operators is very well- 
established [ll| . The feature of FDOs which makes them 
relevant to the context of scale-free transport is their ac- 
tion in Fourier space: 

Fl-ooD^f] = (-ik) a f, F[°°D«f] = (ik) a j, (8) 



where f(k) = J-[f](k). Using property 
the Ricsz FDO satisfies: 



it follows that 



d a f 



d\x\ a 



(k) = -\k\ a f(k). 



(9) 



This result justifies the name "fractional differential op- 
erators" , since they naturally extend the notion of deriva- 
tive to non-integer order. What is most relevant to us is 
that they also provide a suitable generalization of clas- 
sical diffusion. Taking the Fourier transform of Eq. ([5]), 
we obtain: 



dn{k) 



-A a \k\ a n(k), 



(10) 



which reduces to the classical diffusive equation (Eq. ([2])) 
for a = 2 (note that for a = 2 the Riesz FDO yields the 
second derivative operator). This fractional extension 
is relevant for systems in which transport characteristic 
scales diverge. To illustrate this point, let us examine 
the propagator (i.e. the temporal evolution of an ini- 
tial condition n(x,0) = S(x — xo)) of Eq. ([5]). The 
time-derivative of the second moment of the propaga- 
tor is easily related to the transport characteristic length 
scale I. In the case of the diffusive equation (a = 2) the 
propagator is a Gaussian function, 

G(x,t\x ,0) = (47r^)- 1/2 cxp{-(a; -x f/{AAt)} , 

which has a finite second moment that increases linearly 
with time and thus yields a finite I. However, for a < 2, 
the propagator is given by a symmetric Levy distribution 
of order a and behaves like (\x — xq\ » 1): 



G{x,t\x ,0) - C a t- X ' c 



xo\ 



t 1 ^ 



-(l+a) 



(11) 



where C a is a constant [12j. Due to its fat power-law 
tail, the second moment of this propagator is infinite. 
Thus, the transport described by this equation lacks a 
finite characteristic length. This result, together with 
the fact that Levy distributions are stable with res pec t 
to the central limit theorem (as the Gaussian law) [l 21 ] . 
endows FDOs with the physical basis needed to justify 
their use to model scale-free transport. Let us define, for 
convenience, the operator 



a r x — 



-1 



2cos(tt(/3 + l)/2) 



- a D?) . (12) 



Then, the fractional particle flux associated to Eq. © 
can be written as: 



-A ( F c 



Evidently, the operator a F£ satisfies that: 



8 8 a 
— . ( F"- 1 ) = — — 

dx K °° x ' d\xY 



(13) 



(14) 



Eq. (fT5|) is the fractional generalization of the homoge- 
neous Fick law, Eq. ([T]), which is recovered for a = 2 due 
to the fact that ocFf. = d/dx. 

We now come to derive the central result of this Let- 
ter: the inhomogeneous extension of (O that preserves 
GR or, in other words, the scale-free (fractional) general- 
ization of the inhomogeneous Fick law, Eq. ^j. We start 
by giving precise expressions for the GR/LR symmetries 
and then we present a new way to derive the inhomoge- 
neous Fick law (Eq. ([3])). The need for a new approach is 
caused by the failure of the standard methods to derive 
the inhomogeneous Fick law when characteristic length 
scales arc lacking. Finally, we show that the new proce- 
dure works in the fractional case, yielding the fractional 
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generalization of the inhomogencous Fick law. For the 
sake of simplicity, we restrict ourselves to the discussion 
of the Markovian case, but the extension to include non- 
Markovianity is rather straightforward. 

The Continuous-Time Random Walk (CTRW) [H 
provides a suitable framework to derive macroscopic 
transport equations from the microscoscopic dynamics. 
CTRWs are models describing a large number of par- 
ticles whose motion is defined probabilistically using a 
probability density function (pdf). For our purposes, it 
is enough to consider one-dimensional separable mod- 
els defined by two pdfs: a step-size pdf, p(A x ;x), and 
a waiting-time pdf, ip(A t ;x), A t > 0. The joint pdf 
t;(A x ,A t ;x) = p(A x ;x)ip(A t ;x) gives the probability 
that a particle located at x at time t jumps to x + A x at 
time t + A t . As mentioned above, we restrict the discus- 
sion to Markovian problems by selecting a Poisson law for 
the waiting time pdf: ip(A t ;x) = t(x)^ 1 exp (— A t /r(x)), 
where t{x) is the mean waiting time at x. Imposing con- 
servation of the number of particles leads 0, [l5| to the 
Generalized Master Equation (GME) governing the time 
evolution of the density of particles n(x,i): 



dn(x, t) 
dt 



n(x',t) , , , n(x,t) 

——p(x — x ; x )dx -j-^ . (15) 



-(x>) 



t{x) 



The GME is the starting point in the derivation of 
(macroscopic) fluid transport equations, by going to the 
fluid limit in which only long-time, large-distance infor- 
mation is retained. This can be easily done in Fourier 
space by taking the limit of small fc. The Fourier trans- 
form of the GME (|15p with respect to x yields [lj|: 



dh(k, t) 
dt 



n(x',t) 
t(x') 



{p(k- : x')-l)e lkx ' dx' , (16) 



where p(k;x') = p(A; x')e lkA dA is the characteristic 
function of p(A; x'). We introduce now the characteristic 
exponent, A(k;x'), through p(k; x 1 ) = exp A(fc; x'). The 
fluid limit of the GME is obtained then by performing 
the small k approximation p(k] x') ~ 1 + A(fc; .t'), which 
turns Eq. (fT6)) into: 



n ( x '^A(k;x')e ikx 'dx'. (i 7 ) 



dh(k, t) 
dt 



t(x') 



Firstly, let us show how to obtain the homogeneous 
classical (Eq. ^) and fractional (Eq. ([5])) diffusive 
equations from Eq. (TIT]) . We assume that both t(x') 
and A(fc;x') are independent of x' . t becomes then the 
characteristic transport time scale. The characteristic 
length scale is related (if it exists) to the variance <72 
of the step-size pdf, p(A). In the absence of external 
forces, the central limit theorem guides us to choose ei- 
ther a Gaussian or a symmetric Levy distribution for 
p(A) [l^ . The characteristic exponent of the Gaussian 
is A(fc) = — o^fc 2 , which inserted into Eq. (jTTJ) yields the 



(Fourier transform of the) classical diffusive equation (|2|) 
with A = G2/t. The finite variance of the Gaussian is 
related to a finite characteristic length scale I ~ ^fa^. 
In the case of a symmetric Levy pdf, the characteristic 
exponent is [12] A(fc) = -a a \k\ a , with < a < 2. The 
fluid limit of Eq. (JTTJ) then reduces to (the Fourier trans- 
form of) Eq. ([5]). The infinite variance of the Levy pdfs 
implies that / diverges. 

We will now formalize the microscopic symmetries un- 
derlying the inhomogencous Fick law (Eq. ([3])) that we 
discussed briefly at the beginning of this Letter. First, 
we introduce the 1-particle transition rate, T(A; x) = 
p(A; x)/t(x). The LR condition is then expressed as Q 
T(— A, x + A) = T(A, x), V.t. As mentioned, Fick's law 
only requires the less stringent GR condition, 



/OO />OC 
T(— A, x + A)dA = / T(A, x)dA, 
-OO J — OO 



(18) 



as will be shown in the following. 

To derive the inhomogencous Fick law in the classical 
(Gaussian) case, one would typically enforce the sym- 
metry (Eq. (jT5J)) after Taylor expanding the right-hand 
side of Eq. (|15p around x and keeping only the terms 
involving up to the second moment of T . However, a 
Levy pdf of order a < 2 does not have finite moments of 
order equal to or greater than a. Hence, we need to de- 
velop a procedure to find the fluid limit equations while 
imposing GR which does not require an expansion in mo- 
ments. The key point is that, after some straightforward 
manipulations, the symmetry (|18[) can be recast into the 
following, equivalent form: 



/•OO 


rp(fc,x')-ii 


/ — oo 


[ t(x') \ 



c 'dx' = 0, Vfc. 



(19) 



The fluid limit of the symmetry can be expressed in terms 
of the characteristic exponent: 



A(M') ^ 
r(x') 



'dx' = 0, 



(20) 



which is the small fc approximation of Eq. (|19|) . 

It is instructive to rederive the ordinary, inhomoge- 
neous Fick law, Eq. ([3]), using this formalism. Con- 
sider the most general form of the characteristic expo- 
nent of the Gaussian law allowed by the central limit 
theorem 1121: 



A(k;x') = ia(x')k - a 2 (x')k 2 



(21) 



Inserting this expression into (|20[) we immediately find 
that the following relation between a(x') = a(x')/r(x') 
and CT2 (a/) = 02(3/) /t(x') must hold: 



l(fc) = — lfc(T2(fc), 



or, in real space: 



a(x') 



da 2 (x') 
dx' ' 



(22) 



(23) 
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This relation will be recognized as the one invoked ear- 
lier to reduce the full Fokkcr-Planck law to the inhomo- 
geneous Fick law. Using Eqs. (T2Tj) and (f25|) and Fourier 
inverting Eq. (fT7|) . the classical inhomogeneous Fick law 
is recovered: 



dn 



dx 



-0-2{x) 



On 
dx ' 



(24) 



We can use this very same scheme to derive the frac- 
tional version of the inhomogeneous Fick law in the setup 
of Levy pdfs with algebraic tails that are stable according 
to the central limit theorem [T^j • The relevant fractional 
generalization of Eq. (|21[) is now: 

A(k;x') = ia{x')k- a a {x')\k\ a , ag(0,2), (25) 



with a a (x') > 0. Imposing that Eq. ([23)) satisfy again 
GR (Eq. 0) yields ika{k) - \k\ a a^{k) = 0, where 
a(x') = (i(x')/t(x'), a a (x') — o" q (x')/t(x'). The frac- 
tional version of Eq. (l2"2l thus becomes: 



(26) 

{ik) a + (-ik) a we 



a{k) = ±L.fc{k), 

IK 

Using the identity 2cos(7ra/2)|fc| Q = 
can rewrite Eq. (|26|) as 



a(k) = -— a a (k , 

2cos(7ra/2) 



(27) 



and making use of the properties of FDOs under Fourier 
transforms discussed previously: 



a(x') 



(28) 



where the operator ooF 1 " -1 was defined in Eq. (fT2"| . 
Eq. (|28|) is the fractional generalization of the classi- 
cal relation ((23)) . Using this result and Eq. (f25|) in 
Eq. p7[) . and Fourier inverting, we obtain the sought- 
for fractional generalization of the inhomogeneous Fick 
law, dtfi = —d x Tp- a , where the fractional Fick flux is 
now given by: 



[ooF?- 1 (a a n) - (ooF^" 1 ^) n] 



(29) 



Eq. (f2"9"| is the central result of this Letter. It provides 
the keystone on which any description of transport in 
systems lacking characteristic length scales but still sat- 
isfying GR (or LR) should be built. As we mentioned, 
an example could be provided by the propagation of for- 
est fires across an inhomogeneous distribution of trees. 
Another possible application might be the modeling of 
some aspects of radial particle transport in turbulent fu- 
sion plasmas, such as those confined in certain tokamak 



regimes Additional terms might be needed in 



this case to account for additional effects associated to 
radially varying temperatures and external forcing. Note 



also that Eq. (|29|) is very different from the extension 
that one would expect from Eq. (fT3j) in the spirit of Eq. 
©. Namely, 



r_F ;Q ^ -a a (x) 



(30) 



The reason is that fractional derivatives do not satisfy 
the Leibniz rule for the derivative of the product of two 
functions [llj]. Hence, Eq. (|30[) is an equality only in the 
special case a = 2. Needless to say, other fractional gen- 
eralizations (including the one where the flux is defined 
by the right-hand side of Eq. ([30]) ) might also be valid 
under different assumptions. 

Finally, we would like to point out that the scheme 
laid out in this Letter also provides the basis for address- 
ing the derivation of the fractional Fick law in three- 
dimensional inhomogeneous systems. The CTRW con- 
struction can be trivially extended to any dimension. 
However, the formulation in terms of fractional derivative 
operators, even in cases in which rotational invariance re- 
duces the problem to one effective dimension, requires to 
tackle a number of non-trivial technical details, which we 
plan to study in a future work. 
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